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We present results from an analysis of B°{B°) —> p~ using 232 million T(4S') ^ BB decays 
collected with the BABAR detector at the PEP-II asymmetric-energy B Factory at SLAG. We mea¬ 
sure the longitudinal polarization fraction fr ~ 0.978±0.014(stat) Iq 029 (®y®^) ®'iid the GP-violating 
parameters Sl = —0.33 ± 0.24(stat)('(°'°®(syst) and Cl = —0.03 ± 0.18(stat) ±0.09(syst). Using an 
isospin analysis ol B —> pp decays we determine the unitarity triangle parameter a. The solution 
compatible with the Standard Model is a = (100 ± 13)°. 

PACS numbers: 13.25.Hw, 12.15.Hh, ll.30.Er 


In the Standard Model, CP-violating effects in the B- 
meson system arise from a single phase in the Cabibbo- 
Kobayashi-Maskawa (CKM) quark-mixing matrix P|. 
Interference between direct decay and decay after 
mixing in B^^B'^) P~ results in a time-dependent 

decay-rate asymmetry that is sensitive to the angle a = 
arg [—in the unitarity triangle of the CKM 
matrix . This decay proceeds mainly through a 6 —> uud 
tree diagram. The presence of penguin loop contributions 
introduces additional phases that shift the experimen¬ 
tally measurable parameter Oeff away from the value of 
a. However, measurements of the B'^ P^P^ branching 

fraction and the upper limit for B^ —> p°p° Si show 
that the penguin contribution in B pp is small with 
respect to the leading tree diagram, and Sapp = Oeff —a is 
constrained at ±11° at Icr i. This Letter presents an up¬ 
date of the time-dependent analysis of B^{B^) p^ p~ 

and measurement of the CKM angle a reported in i . 

The CP analysis of B decays to p'^ p~ is complicated by 
the presence of a mode with longitudinal polarization and 
two with transverse polarizations. The longitudinal mode 
is CP even, while the transverse modes contain CP-even 
and CP-odd states. Empirically, the decay is observed to 
be dominated by the longitudinal polarization 0, with 
a fraction Jl defined by the fraction of the helicity zero 
state in the decay. The angular distribution is 

_ fL _= ( 1 ) 

Tdcos 9id cos 62 

^ /l cos^ 9 i cos^ ^2 ± ^ (1 — /l) sin^ 9 i sin^ 62 

where 0i=i,2 is the angle between the tt^ momentum and 
the direction opposite the B^ in the p rest frame, and 
we have integrated over the angle between the p decay 
planes. 

The analysis reported here is improved over our earlier 
publication 0 by a change in selection requirements re¬ 
sulting in an increased signal efficiency; introduction of 
a signal time dependence that accounts for possible mis- 
reconstruction; and use of a more detailed background 
model. This measurement uses 232 million 7’(4iS') ^ BB 
decays collected with the BABAR detector at the PEP- 
II asymmetric-energy B Factory at SLAG. 


We reconstruct B^{B°) P^P~ candidates (Prec) 
from combinations of two charged tracks and two tt^ can¬ 
didates. We require that both tracks have particle identi¬ 
fication information inconsistent with the electron, kaon, 
and proton hypotheses. The 7r° candidates are formed 
from pairs of photons each of which has a measured en¬ 
ergy greater than 50 MeV. The reconstructed 7r° mass 
must satisfy 0.10 < m-yj < 0.16 GeV/c^. The mass of the 
p candidates must satisfy 0.5 < m^±^o < 1.0 GeV/c^. 
When multiple B candidates can be formed, we select 
the one that minimizes the sum of — m^o)'^ where 
iriTrO is the true tt^ mass. If more than one candidate has 
the same 7r° mesons, we select one at random. 

Gombinatorial backgrounds dominate near | cos 9i \ = 
1, and backgrounds from B decays tend to concentrate 
at negative values of cos 9i. We reduce these backgrounds 
with the requirement —0.90 < cosOi < 0.98. 

Gontinuum qq {q = u,d,s,c) events are the 

dominant background. This background is reduced by 
requiring that IcosPr^l < 0.8, where Btr is the an¬ 
gle between the B thrust axis and that of the rest of 
the event, ROE. The thrust axis of the B is the direc¬ 
tion which maximizes the longitudinal momenta of the 
particles in the B candidate. To distinguish signal from 
continuum we use a neural network (A/") to combine ten 
discriminating variables: the event shape variables that 
are used in the Fisher discriminant in Ref [^; the cosine 
of the angle between the direction of the B and the col¬ 
lision axis (z) in the e+e” center-of-mass (GM) frame; 
the cosine of the angle between the B thrust axis and 
the z axis, |cosPt/j|; the decay angle of each 7r° (de¬ 
fined in analogy to the p decay angle, 9i); and the sum 
of transverse momenta in the ROE relative to the z axis. 

Signal events are identified kinematically using two 
variables, the difference AE between the GM energy 
of the B candidate and \/sj2, and the beam-energy- 
substituted mass rriES = \/ (s/2 ± • psY/Ef — p^, 

where ^/s is the total GM energy. The B momen¬ 
tum pb and four-momentum of the initial state (Pi,Pi) 
are defined in the laboratory frame. We accept can¬ 
didates that satisfy 5.23 < toes < 5-29 GeV/c^ and 
—0.12 < AE < 0.15 GeV. The asymmetric AE selection 
reduces background from higher-multiplicity B decays. 
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To study the time-dependent asymmetry one needs to 
measure the proper-time difference, At, between the two 
B decays in the event, and to determine the flavor of 
the other B meson (Btag)- We calculate At from the 
measured separation Az between the Srec and Stag decay 
vertices 0- We determine the i?rec vertex from the two 
charged-pion tracks in its decay. The Stag decay vertex 
is obtained by fitting the other tracks in the event, with 
constraints from the Brec momentum and the beam-spot 
location. The RMS resolution on At is 1.1 ps. We only 
use events that satisfy |At| < 20 ps and for which the 
error on At less than 2.5 ps. The flavor of the i?tag meson 
is determined with a multivariate technique 0 that has 
a total effective tagging efficiency of (29.9 ± 0.5)%. 

Signal candidates may pass the selection requirement 
even if one or more of the pions assigned to the p'^ p~ 
state belongs to the other B in the event. These self- 
cross-feed (SCF) candidates constitute 50% (26%) of the 
accepted signal for /^ = 1 (/^ = 0). The majority of 
SCF events have both charged pions from the p'^p~ final 
state, and unbiased CP information (correct-track SCF). 
There is a SCF component (14% of the signal) where at 
least one track in Spec is from the rest of the event. These 
wrong track events have biased CP information, and are 
treated separately for the CP result. The probability 
density function (PDF) describing wrong track events is 
used only in determining the signal yield and polariza¬ 
tion. A systematic error is assigned to the CP results 
from this type of signal event. 

We obtain a sample of 68703 events that enter a 
maximum-likelihood fit. These events are dominated by 
backgrounds: roughly 92% from qq and 7% from BB 
events. The remaining 1% of events is signal. We dis¬ 
tinguish the following candidate types: (i) correctly re¬ 
constructed signal; (ii) SCF signal, split into correct and 
wrong track parts; (iii) charm B^ background {b ^ c); 
(iv) charm background {b c); (v) charmless B 
backgrounds; and (vi) continuum background. The dom¬ 
inant charmless backgrounds are B decays to pn, (citt)^, 
(ai7r)°, and longitudinally polarized aip final states. For 
these decays we use the inclusive branching fractions 
(in units of 10-6), 34 ± 4 0, 42 ± 42, 42 ± 6 0 and 
100 ± 100, respectively. The corresponding expected 
number of events in the sample are 82 ± 13, 87 ± 87, 
65 ± 9, and 202 ± 202. We also account for contribu¬ 
tions from higher kaon resonances (112 ± 112 events) and 
p'^p^ (82 ± 19 events). In addition we expect 2551 ± 510 
(1316 ± 263) charged (neutral) B decays to final states 
containing charm mesons. The R-background decays are 
included as separate components in the fit. 

Each candidate is described with the eight R^ec kine¬ 
matic variables: toes, AE, the m,r±- 7 r 0 and cos 9i values 
of the two p mesons, At, and JV. For each fit component, 
we construct a PDF that is the product of PDFs for 
these variables, neglecting correlations. This introduces 
a fit bias that is corrected with the use of Monte Carlo 


(MC) simulation. The continuum-background yield and 
its PDF parameters for toes, AE, cos9i, and A/” are 
floated in the fit to data. The continuum m^±^o dis¬ 
tribution is described by a Breit-Wigner and polynomial 
shape, and is derived from toes and AE data sidebands. 
For all other fit components the PDFs are extracted from 
high-statistics MC samples. The cos^i distributions for 
the background are described by a non-parametric (NP) 
PDF derived from the MC samples, as the detector ac¬ 
ceptance and selection modify the known vector-meson 
decay distribution. The true signal distribution is given 
by Eq. ^multiplied by an acceptance function determined 
from signal MC samples, whereas SCE signal is modeled 
using NP PDFs. 

The signal decay-rate distribution for both polariza¬ 
tions /+(/-) for Rtag= B'^ (R6) is given by 

g-|At|/T 

f±{At) = —--[1 ± S'sin(ATOdAf) =F C'cos(ATOdAf)], 

where r is the mean B^ lifetime, Ato^ is the B^B^ mix¬ 
ing frequency, and S = Sl or St and C = Cl or Ct 
are the CP-asymmetry parameters for the longitudinally 
and transversely polarized signal. The parameters S and 
C describe P-mixing induced and direct CP violation, 
respectively. S and C for the longitudinally polarized 
wrong-track signal are fixed to zero. The At PDF takes 
into account incorrect tags and is convolved with the res¬ 
olution function described below. Since /l is approxi¬ 
mately 1, the fit has no sensitivity to either St or Ct- 
We set these parameters to zero and vary them in the 
evaluation of systematic uncertainties. 

The signal At resolution function consists of three 
Gaussians (^90% core, ^9% tail, ^1% outliers), and 
takes into account the per-event error on At from the 
vertex fit. The resolution is parameterized using a large 
sample of fully reconstructed hadronic B decays 0 . For 
wrong-track SCF we replace the P-meson lifetime by an 
effective lifetime obtained from MC simulation to account 
for the difference in the resolution. The nominal At dis¬ 
tribution for the P backgrounds is a NP representation 
of the MC samples; in the study of systematic errors we 
replace this model with the one used for signal. The reso¬ 
lution for continuum background is described by the sum 
of three Gaussian distributions whose parameters are de¬ 
termined from data. 

We perform an unbinned extended maximum likeli¬ 
hood fit. The results of the fit are 617 ± 52 signal 
events, after correction of a 68 event fit bias, with Jl = 
0.978±0.014, Sl = -0.33±0.24 and Cl = -0.03±0.18. 
The measured signal yield, polarization, and CP param¬ 
eters are in agreement with our earlier publication 0 , 
with significantly improved precision. Figure ^ shows 
distributions of toes, AP, cos 9i and TO,r±7rO for the high¬ 
est purity tagged events with a loose requirement on 
J\f. The plot of Toes contains 14% of the signal and 
1.5% of the background. For the other plots there is an 
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FIG. 1: The distributions for the highest purity tagged events 
for the variables mss (a), AU (b), cosine of the p helicity 
angle (c) and m^±^o(d). The dotted lines are the sum of 
backgrounds and the solid lines are the full PDF. 



FIG. 2: The At distribution for a sample of events enriched in 
signal for 5° (a) and (b) tagged events. The dotted lines 
are the sum of backgrounds and the solid lines are the sum of 
signal and backgrounds. The time-dependent CP asymmetry 
(see text) is shown in (c), where the curve is the measured 
asymmetry. 


added constraint that tties > 5.27 GeV/c^; these require¬ 
ments retain 11.5% of the signal and 0.4% of the back¬ 
ground. Figure 121 shows the At distribution for and 
tagged events. The time-dependent decay-rate asym¬ 
metry [A(At) — A(At)]/[A(At) + A(At)] is also shown, 
where N (N) is the decay-rate for B° {B°) tagged events. 

We have studied possible sources of systematic uncer¬ 
tainties on /l, Sl and Cl- The dominant uncertain¬ 


ties for /l come from floating the B background yields 
(Ig;™), non-resonant events (0.015) and fit bias (0.01). 
The dominant systematic uncertainty on the CP results 
comes from the uncertainty in the i?-background branch¬ 
ing ratios. This results in a shift on Sl (Cl), as large as 
-o ’12 ( -o’oos)- Additional uncertainties on the CP results 
come from possible CP violation in the B background, 
calculated as in Ref. 0 . We allow for a CP asymmetry up 
to 20% in B decays to final states with charm, resulting 
in an uncertainty of 0.027 (0.045) on Sl {Cl)- Allowing 
for possible CP violation in the transverse polarization 
results in an uncertainty of 0.02 (lo;o?g) on Sl {Cl)- 
We estimate the systematic error on our CP results from 
neglecting the interference between B°{B^) p'^p~ and 

other 47r final states: B aiir, pinr^ and B 7r7r7r°7r°. 
Strong phases and CP content of the interfering states are 
varied between zero and maximum using uniform prior 
distributions, and the RMS deviation of the parameters 
from nominal is taken as the systematic error; this is 
found to be 0.02 on Sl and Cl- Other contributions 
that are large include knowledge of the vertex detector 
alignment 0.034 (0.005) on Sl {Cl), and possible CP vi¬ 
olation in the doubly-Cabibbo-suppressed decays on the 
tag side of the event M- We allow CP violation in the 
wrong-track SCF to vary between —1 and -|-1, which re¬ 
sults in changes of 0.007 (0.012) in Sl {Cl)- The nominal 
fit does not account for non-resonant background. If we 
add a non-resonant component oi B —> events to 

the likelihood, we fit 83 ± 59 non-resonant events and 
observe only a (6 ± 4)% drop in signal yield. This effect 
is included in our total systematic uncertainty. Possible 
contributions from cr(400)7r°7r° decays are neglected due 
to the small reconstruction efficiency (0.4%). Our results 
are 

h = 0.978 ±0.014(stat) 1^29 (syst), 

Sl = —0.33± 0.24(stat)^Q°4(syst), 

Cl = -0.03±0.18(stat)±0.09(syst), 

where the correlation between Sl and Cl is —0.042. 

We constrain the CKM angle a from an isospin analysis 
of B ^ pp. The inputs to the isospin analysis are 
the amplitudes of the CP-even longitudinal polarization 
of the pp final state, as well as the measured values of Sl 
and Cl for P°(P°) —> p'^p~ - We use the measurements 
of /l, Sl and Cl presented here; the branching fraction 
of B^ p'^ p~ from 0 , which uses information from 
0; the combined branching fraction and /i for B —> 
p+p° from Ref. 0; the central value corresponding to 
the upper limit oiB{B p°p°) from Ref. 0 . We ignore 
electroweak penguins and possible 7 = 1 amplitudes . 

To interpret our results in terms of a constraint on 
a from the isospin relations, we construct a that in¬ 
cludes the measured quantities expressed as the lengths 
of the sides of the isospin triangles and we determine 
the minimum Xo- As the isospin triangles do not close 
with the current central values of the branching ratios. 
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we have adopted a toy MC techniques to compute the 
conhdence level (CL) on a; our method is similar to the 
approach proposed in Ref. 0 . For each value of a, 
scanned between 0 and 180°, we determine the difference 
Axqata(®) between the minimum of x^(a) and Xo- We 
then generate MC experiments around the central val¬ 
ues obtained from the fit to data with the given value 
of a and we apply the same procedure. The fraction 
of these experiments in which Ax|j( 2 (a) is smaller than 
Axqata(®) interpreted as the CL on a. FigureOlshows 
1 — CL for a obtained from this method. Selecting the 
solution closest to the CKM combined fit average 
we find a = 100° ± 13°, where the error is dominated by 
6 app which is ±11° at lu. The 90% CL allowed interval 
for a is between 79° and 123°. 



FIG. 3: CL on a obtained from the isospin analysis with 
the statistical method described in M- The dashed lines 
correspond to the 68% (top) and 90% (bottom) CL intervals. 


In summary we have improved the measurement of the 
CP-violating parameters Sl and Cl in P°(P°) ^ 
using a data-sample 2.6 times larger than that in Ref. |j|. 
We do not observe mixing-induced or direct CP viola¬ 
tion. We derive a model-independent measurement of 
the CKM angle a, which is the most precise to date. 
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